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ABSTRACT 

We describe a model-independent integral constraint which defines an upper limit to the allowed degree of 
evolution in the ratio of black hole (BH) mass to host galaxy luminosity or mass, as a function of redshift. 
Essentially, if the BH/host ratio is excessive at redshift z, then it would imply that the total mass density in BHs 
above some Mg^j 1 is larger at that redshift than at z = 0, which is impossible. This argument requires no knowl- 
edge of host or BH properties, only a lower limit to the observed luminosity density in the brightest galaxies at 
some z- We calculate this upper limit from a variety of luminosity and mass functions in different bands from 
redshifts z = 0-2. We show that it is consistent with passive evolution of spheroid populations (with a fixed 
M BH /M host relation) in all cases, and provides tighter constraints than have generally been obtained previously, 
ruling out at > 6a observational and theoretical estimates suggesting thatMBH/Mh os t was significantly larger at 
high redshifts than locally, although relatively weak (factor < 2 by z = 2) evolution is still allowed. We discuss 
a variety of possible "loopholes" or changes in the BH/host populations and correlations, and show that they 
typically lower the upper limits and strengthen our conclusions. 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 
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1. INTRODUCTION 

Recent discoveries of tight correlations between the 
masses of supermassive black holes (BHs) in the cen- 
ters of nearby galaxie s and either the l uminosity (e.g., 



mass (Magorrian et al. 



Kormendv & Richstone 

19981). or velocity dispersion jFerrarese & Merrittl 120001: 



Gebha rdt et al J 12000) of their host spheroids demonstrate a 
fundamental link between the growth of BHs and galaxy 
formation. Determining the evolution of these correlations 
with redshift is cri tical for informing analytical model s (e.g. , 
Silk &Reesi 119981) an d simulations (Pi M atteo et all l2005t 
Robertson et al . 2006) which follow the co-formation of BHs 



and bulges, as well as theories which relate the evolution and 
statistics of BH formation and quasar activity to galaxy merg- 
ers (e.g., Hopkins et al. 2006a d) and to the remnant spheroid 
population (Hopkins et al. 2006b). Likewise, the significance 
of observa tions tracing t he buildup of spheroid populations 
(e.g..lCowie etalJfl996t IWolf et alJl2003t iBell et alJ l2004a: 
Bund v et alJl2005bl) and associations betwee n spheroids in 
formation, mergers, and quasar hosts (e.g., Hopki ns et alJ 
l2006d) depends on understanding the evolution of BH/host 
correlations. 

Efforts to directly infer these correlations at redshifts z > 
are, however, difficult and limited by the small num- 
bers of objects with observable hosts. Consequently, dif- 
ferent groups have reached seemingly contradict ory conclu- 
sions, with 0[III] v elocity dispersion JShields et alJ l2003h 
and BH clustering ( Adelberger & Steidell 120051) measure- 
ments finding no evolution in these corr elations out to red- 
shifts z ~ 4, while C O velocity dispersion ( Walter et all2004t 
Shield set alJl20 05) and host 7?-band luminosity (Pe ng et all 
2006) observations imply significant evolution (with BHs 
becoming over-massive relative to their hosts) occurring at 
z > 1 . Spectral template fitting ( iTreu et alJl2004l IWoo et alJ 
2006) suggests strong evolution at much lower redshifts z = 
0.37, and comparison of radio quasar and galaxy popula- 
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tions (McL ureet all 2006) indicates evolution of the form 
oc (1 +z) 2 . In addition to the difficulty of identifying host 
properties at these redshifts, the measu rements generally rely 
upon virial line-width relations (e.g., Kaspi et al. 20051 an d 
references therein) to determine BH masses, which have not 
been tested at high redshifts and require that the BH be visible 
as an optical, broad-line AGN. 

It is therefore of interest to develop additional constraints on 
the evolution of BH/host correlations, especially ones which 
are independent of systematics in host measurements and se- 
lection effeetsJri_considering only AGN populations. For ex- 
ample. lMerloni et aflll2004l) compare the integrated quasar lu- 
minosity and stellar mass densities, and infer a much weaker 
evolution of the form oc (1 +z) 5 , but this is subject to sys- 
tematic uncertainties in quasar radiative efficiencies and (as 
discussed below) in observational completeness and galaxy 
type segregation at low luminosities. In this paper, we formu- 
late and apply a model-independent integral constraint which 
provides an upper limit to evolution in these correlations as 
a function of redshift, and show that it rules out various pro- 
posed forms of evolution at high significance. We adopt a 
fi M = 0.3, Cl\ = 0.7, H = 70kms _1 Mpc" 1 cosmology, but the 
choice of Hq factors out in our analysis and the limits depend 
only weakly on Qm and J7a- 

2. THE INTEGRAL CONSTRAINT 

If there exists some mean relationship between BH mass 
(Mbh) and host mass or observed luminosity in a (rest-frame) 
band (Lh), we can parameterize it with the coefficient /z(z), 

(M BH )=M(z)^- (1) 
For simplicity in what follows, we adopt a linear correlation, 
but the argument can be generalized trivially for any func- 
tional dependence Mbh = A i (z)/z(^h)- The term Lh can rep- 
resent either total or strictly bulge/spheroid luminosity. Inte- 
grating the space density of both sides of Equation Q down 



to some minimum = /i(z)L™ 



i.e. 



M BH dn(M BH ) = n(z) 

M%2 Jli 



L h dn(L h ) 



(2) 
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Hopkins et al. 



gives an equivalent relation between BH mass density and 
galaxy (or spheroid) luminosity density (jh), 

Pbh(M bh > MZ , z) = M(z) h(L h > MZ/fiiz) , z) , (3) 



MEff/M.) 
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with the z = equivalent Pbh(M B h > Mgn , 0) = /i(O) ;h(L h > 
M££?//i(0) , 0). Note that scatter in the M BH -L h relationship 
will introduce an additional multiplicative factor here, but this 
will generally cancel in our analysis (discussed further be- 
low). We divide the relations to obtain 

Pbh(M bh > M™ n , z) fi(z) MU > Mgfr/Kz) , z) 



Pbh(M bh > M™ 1 , 0) 



M(0); h (L h >M™ n M0),0) 
ML h >L /T(z),z) 



= r( Z ) 
<i. 



j h (L h > L , 0) 



(4) 



where the second equality comes from substitution with the 
definitions T(z) = m(z)/m(0) a nd L = M™ n //x(0 ). Note that 
this notation for T(z) is used bv lPeng et al.l (l2C)06). but in their 
case refers to individual systems, rather than the evolution of 
the mean relation. The third equality comes from the physical 
fact that it is not possible to destroy BH mass - the total mass 
above a given Mg^ can never decrease from z to the present. 
This is true for every M'^ and corresponding Lo, which im- 
plies an upper limit to T(z) for arbitrary Lo, 



, ;'h(L h > L /r max (z) , z) _ 

X w . . —I, 

j h (L h > Lo , z = 0) 



(5) 



where again the jh are defined in the same rest-frame band. 
Since j'h(Lh > x) is by definition a monotonic decreasing func- 
tion of x, the left hand side of Equation Q is necessarily a 
monotonic increasing function of T(z), and the upper limit 
for a given Lo is well-defined. Any other functional depen- 
dence where BH mass increases with host luminosity or mass 
at fixed redshift will give a qualitatively identical constraint 
(it just will not lend itself to such convenient notation). 

Given a luminosity function (LF) or luminosity density 
above some luminosity jh(Lh > Lo) at z = and some higher 
redshift, we can determine the upper limit r max (z) for ev- 
ery minimum luminosity Lo. Figure illustrates how this 
depends on the "dummy" variable Lo, with an example of 
r m ax(z, Lo) from Equation 0. This is heuristic for now, but 
in deta il we use the ^ -band observed local spheroid/red galaxy 
LF o flBell et aflj2003l) to calculate j h (L h > L , z = 0), and the 
z = 0.2-0.4 and z = 1.0-1.2CO MBO-17 fi-band spheroid/red 
galaxy LFs from |Faberetal] d2005b to calculate ih(Lu > 
L /r(z) , z) (see also lBell et al\m^ - $m\mer et al.l (2005) 
for details and the nearly negligible B-g correction). We cal- 
culate the upper limit defined by Equation[3](r max ) for every 
minimum luminosity L ( > at which the binned LF data exists 
(squares), with a \a upper limit (vertical error bar) estimated 
from the observational uncertainty in jh(Lh > Lq /T(z)) using 
Equation Q to determine ST m:ix /5jh (technically we also add 
the errors in jh from z = in quadrature, but this is in all cases 
dominated by the jh(z > 0) errors). 

How strong the constraint is depends on which Lo we inte- 
grate down to. For example, integrating to low Lo effectively 
compares the cumulative BH mass densities at two redshifts, 
whereas integrating to a relatively high Lo compares the BH 
mass density in only the most massive BHs. Both must be 
lower at z than locally, but if e.g. the majority of the most 
massive BHs are in place by z ~ 1 (such that their mass den- 
sity is already comparable to that observed locally) while only 
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FIG. 1. — Upper limit to the allowed evolution in the proportionality be- 
tween BH mass and host (spheroid or total) B-band luminosity (fj,(z) / /4z = 0), 
where Mbh = as a function of the (arbitrary) minimum luminosity 

(Lo) to which we integrate the luminosity density in Equation at z = 0.3 
and z = 1. 1 from comparison o f the observed luminosity density in bright 
spheroids I Bell et al. 2003; Faber et al. 2005). Points show the upper limits, 
eiTor bars the uncertainty (dominated by statistical uncertainty in the total 
luminosity density observed at L > Lq). For comparison, the corresponding 
Mg^j 1 = /i(0)Lo is also shown, although the actual values of fi(0) and Lq are 
not important. If T(z) = fi(z)/fJ.(z = 0) were greater than the limit shown at 
any Lq, then the BH mass density above the corresponding Mg 1 ^ 11 would be 
larger at z than it is locally. Equation |4j must hold for all Lo, so the lowest 
statistically significant observational limit is the appropriate constraint. 



a small fraction of low mass BHs exist (such that the total BH 
mass density may be substantially lower than that at z = 0), 
then the former integration will give a much stricter limit 
r max . But since it is an upper limit which must hold above 
all minimum masses or luminosities Lo, the lowest value of 
Lmax is the proper constraint. 

This tends to derive from the highest luminosities because 
observations indicate that the stellar mass of the largest ob- 
jects is largely for med by z ~ 1 ("cosmic downsizing," e.g., 
ICowie et alJl!996l) . meaning that the luminosity density in 
them is constant or rises with z, as the stellar populations are 
younger and brighter at higher redshift. As long as the num- 
ber of objects at the highest masses does not fall dramatically 
with z, this trend in the constraints on T(z) is preserved. In 
other words, if a significant fraction of the most massive ob- 
jects are formed by z ~ 1 and their rest-frame luminosity is 
typically brighter at that time, then the BH mass to host lu- 
minosity ratio cannot also increase, or the mass density in the 
most massive BHs would be much higher at z ~ 1 than it is 
today. 

The observations are more complete and are dominated by 
spheroids at these luminosities, which reduces our uncertain- 
ties and also means that it does not matter whether the appro- 
priate correlation is with spheroid or total galaxy light/mass. 
However, at the brightest luminosities, the number of ob- 
served objects is small and statistical uncertainty is larger. 
Therefore, in what follows, we take as our constraint the value 
of r max which has the minimum la upper limit. This is a con- 
servative assumption, as higher Lo will give lower (but less 
well-constrained) limits and the cumulative significance of the 
upper limit (i.e. cumulative la upper error with respect to the 
entire data set) will always be smaller than that we show. 

With a well-defined r max from any observed sample at a 
given redshift, we now consider r max (z). Figure|2]shows this 
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upper limit as a function of redshift, from a large number of 
different observed LFs, in B, r, and K bands as well as es- 
timated directly from stellar mass functions. We calculate 
r m ax (squares) and its error as above in each case, and also 
plot (horizontal error bars) the redshift range over which each 
binned LF is measured. Because it appears locally that the ap- 
propriate correlation is between Mbh and the bulge/spheroid 
luminosity or mass, we consider only samples which have 
measured LFs for both early and late types, and derive our 
constraints only from the early type luminosity density com- 
parison. We discuss this below, however, and find it makes 
no difference to our conclusions. For each band, we compare 
to the best-determined local luminosity density as a function 
of the minimum luminosity Lq. The resulting constraints are 
generally strongest from optical LFs, because these typically 
involve the largest samples and thus have the best statistics at 
the bright end. The mass-function constraints, by compari- 
son, are similar where the samples are large enough to probe 
comparable number densities, but generally derive from much 
smaller samples. 

For comparison, we plot r max (z) expected from passive 
evolution, assuming a constant BH to host stellar mass ra- 
tio (Mbh/M*) with the spheroid luminosity and colors pas- 
sively evolving from a formation redshift Zf = 1, 2, 3, 5 (long 
dashed, dotted, short dashed, solid lines, respectively) at solar 
metall icity according to the stella r population synthesis mod- 
els of Bruzual & Chariot (2003). We also plot (open black 
circles) the r ecent estimates o f the evolution in these corre- 
lations from Peng et al. (2006). The authors directly mea- 
sure the correlation with R; we use their adopted model to 
convert to mass, B, and K band estimates. Note that they 
assume passive evolution with zj = 5; as they point out, al- 
lowing for star formation or younger populations will push 
their estimates of T(z) to higher values. In the stellar mass 
panel, we also sho w (open blac k points) the observational 
estimates of Shields e\ alj (|2Q0 3[ cro sses'). IWoo ^yilJ (2006. 
square), Adelberger & Steidel (2005, star), and McLure et al. 
( 2006, diamonds). The dotted line in each panel shows Mbh oc 
(1 +z) 2 M* with Zf = 5, a rough upper envelope of these obser- 
vations as they would manifest in other bands. 

3. DISCUSSION 

Our model-independent integral measure provides a con- 
straint on the evolution in the correlation between BH mass 
and host luminosity or mass. The upper limits we derive are 
competitive with or stronger than those from direct, but dif- 
ficult observations of high-redshift BHs and their hosts, and 
can be applied in any band (including those for which direct 
estimates of these correlations at z > do not exist). Our 
upper limits are consistent with passive evolution (with rea- 
sonable formation redshifts z/ < 6) at all redshifts and in ev- 
ery band, as well as hydrodynamic simulat ions of BH and 
spher oid co-formation at different redshifts dRobertson et alJ 
2006), and vari ous estimates of these correlations at low 
redshifts (z < 1. IShields et alJl2003t IPeiTg et al]l2006l) . The 
weak evolution inferred by Merl oni et alJ (|2004) cannot be 
ruled out for sufficiently low z/ < 3, but the mass func- 
tion data prefer a no-evolution case at ~ 2 a. However, the 
best-fit evolution e stimates (open points in Figure E) by e.g . 
Treu et al.1 J200l: IShieldsetaD J2005h : iPeng et all (2006); 
McLure et al.1 (120061) i lWoo et all (120061) are ruled out at > 6a 
with respect to the upper limits we measure over each cor- 
responding redshift interval. (Note that this does not neces- 
sarily mean their measurements are inconsistent with our up- 



per limits at 6a, as the error bars in the direct observations 
are large and generally inconsistent with passive evolution at 
only ~ 2-3cr.) It should be noted, however, that there is a 
factor ~ 2 systematic normalization uncertainty in the virial 
BH mass est imators the s e auth ors adopt, and lowering e.g. the 
estimates of Peng et al. (2006) by such a factor (i.e. allowing 
for a factor of ~ 2 increase in Mbh/M* by z ~ 2) makes them 
consistent with our upper limits. 

Given this apparent contradiction, it is natural to ask 
whether there are effects which could increase the upper lim- 
its (r max ) in Figure|2] Observational incompleteness will only 
lead to a higher "true" luminosity density at some z, and lower 
r ma x (assuming the luminosity density at z = is not substan- 
tially less complete than that at higher redshift). If we have 
somehow severely overestimated the mass in large spheroids 
at high z, i.e. if these red systems are really all dusty starbursts, 
this could increase r max , but many of the samples shown in 
Figure |3 are morphologically, not color, selected, and obser- 
vations find that only a small fraction < 10% of the luminosity 
density in red galaxies at z ~ 1 comes from dusty, in trinsically 
blue so urces (too small to change our conclusions) ( B ell et all 
l2004bl). However, this m ay rise to ~ 30-50% at 1 < z < 2 
(e.g- ICimatti et al.l20 02). shifting t he limits at these reds hifts 
from the color-select ed samples of Fontana et al. (2003' and 
Giallong o et al.l (12005) up by a factor ~ 1 .5 - 2. 

Furthermore, the most meaningful upper limits come from 
high luminosities (typically ~ 0.5 - 1 mag above the LF break 
at each z), where a significant fraction of the most massive 
galaxies are in place by z ~ 1 and this constraint is most sen- 
sitive to a change in /i(z), but also where the observations are 
most well-constrained. For this reason, it makes no difference 
which local LF we use to compare in Equation (0, as they 
all agree we ll at th ese luminosities (see e.g. Figures 8 & 9 
of Bell et al. (2003)) . We have al s o dire c tly tested t h is wit h 
the local LFs in iKochanek et all (120011): iBell et alJ (120031) : 
iBlanton et ail (120031) : iNakamura et al.1 J2003I) and find the re- 
sults are similar. Moreover, it makes no difference whether 
or not we use the cumulative (all galaxy) LFs (i.e. assume a 
correlation between Mbh and total galaxy light or mass) in- 
stead of spheroid LFs as above, since spheroids dominate the 
luminosity density above these bright luminosities. In any 
case, the contribution of late-types tends to evolve in a simi- 
lar fashion to that of early type s (see e.g. the observations in 
Figure 13 and D rorv et al.ll2004l) . In even more detail, allow- 
ing for late types but accounting for their bulges (with more 
detailed mo rphology-de pendent LFs) makes negligible differ- 
ence, as e.g. iKoo et al.l d2005 ) find most massive bulges are 
already in place at z ~ 1 , giving a similar constraint but with 
more objects and thus a stronger la upper limit. 

Changing the form of the correlation also does not al- 
ter our results. Our derivation of Eqn. [5] can be performed 
for any correlation of the form Mbh = fi(.z)f(Lh) or Mbh = 
/[/i(z)LhL and any relationship where Mbh is an increasing 
function of host mass or luminosity will give a similar result. 
Furthermore, all non-linear forms of the correlation /(Lh) 
which have been proposed observationally or theoretically 
steepen at high luminosity (high Mbh), increasing the im- 
plied mass of the most massive BHs and their relative weight 
in the BH mass density, and directly considerin g these non- 
linear (iKormendv & Richstondll995tlMarconi & Huntll2003t 
Haring & Rix 2004) or even log-quadratic (Wvithe 2006) cor- 
relations in every case lowers r max . 

But what if the shape or slope of /(Lh) changes with red- 
shift? This would be particularly interesting, in the implica- 
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FIG. 2. — Upper limit to the allowed evolution in the proportionality between BH and host (spheroid or total) luminosity (Mbh = Mz)^h) as a function of 
redshift, based on t he observed lumin osity density or stellar mass densit y as in F igure[Tl(sauares, with upper lcr error from ob servational uncertainti es). Limits 
in fl-band are fromlBell et alj (SDSS+2M ASS, local zero-point). [Blanton etal . 1 2003) (SDSS, green). limetaTI E002 ) (DEEP GS S, oran ge), IBell et alj 

i 2004a1: lWiHmer et alJ<2005ft:IFaber et alj|2 005) (COMBO-17, blue, D EEP2, red),IDahlen et alJ<2T)05ft (GOO DS, purple). Giall ongo etaflf 2005) (COM BO-17 , 
cyan). In r-band from Nakamur a et alJ 12003ft ("SPSS, l ocal zero-p oint), Blanton et al. 1 20031) f green). IZucca etalj J2005) (VVDS, blue), Gabasch e^uj j^tOffi 
(EPF cyan), IChen et alj 12003ft ( LCIR Survey, orange) , IWolfeFal] 120031) (COMBO-17^ red). In g-band from lBell et alTl2003T) (loca l zero-p oint), ICowie et all 
119961) ( Hawaii Deep Fields, red), Pozzetti et al. 1 2003) (K20, blue), Dahlen et al. 1 2005) (purple). Mass function s are from Bell et al. 1 2003$ (local zero-point), 
IBundv et all 12005a b) (GOODS+DEEP2, red), |Franceschini et al. i2006D (IRAC/GOODS, cyan). IFontana et alj (K20, blue). All bands are rest-frame. 

Curves show the expectation of passive evolution tracks (with constant Mbh A/*), with fo rmation redsh ifts zf = 1, 2, 3, 5 (long-dashed, dot-dashed, s hort-da shed, 
solid, respectively). Open black points show the recent estimates of Shields et al. 1 2003, crosses). Woo et al. 1 2006, square), Adelberger & Steidel 1 2005, star), 
IPeng et all |2*006. circles), McLure et al. 1 2006, diamonds); dotted line shows Mbh <x (1 +z) 2 M* with z/ = 5. The upper limits are consistent with passive 
evolution in all cases, but strongly rule out a substantially larger Mbh/M* at high redshift. 



tion that the ratio Mbh/M* evolves differently for systems of 
different mass. The exact effect on our derivation depends on 
the "pivot point" about which the relationship slope changes. 
For example, given that most of the Mbh _ measurements 
in Figure |2 have been restricted to high-M BH , keeping these 
points fixed but finding a steeper slope at high redshift would 
imply a lower observed T(z) at lower masses, which could re- 
lieve tension between our upper limits and the observations, 
although it would also raise the implied T(z) at higher masses. 
In such a case, our integral constraint still obtains, it simply 
must be applied as a function of Mf^j 1 or Lq, as T(z) becomes 
a function of mass as well as redshift. The limits shown in 
Figure |2] would then refer to the specific MJ^j 1 to which we 
perform our integration (generally ~ 10 8 - 10 9 M Q ). 

Allowing for scatter in the BH/host relationship /(Xh) like- 
wise does not alter our results. Non-zero scatter a{z) will in- 
troduce a term ~ exp{0.5 In 2 10(cr 2 (z)-cr 2 (0))} in the left- 
hand side of Equation Q (the appropriate exact factor will 



be a non-trivial function of Lrj, but since the relevant con- 
straints come from somewhat above the LF break, this factor 
is accurate to ~ 10% and the qualitative behavior is identi- 
cal). A larger scatter implies more objects at the highest Mbh 
for a given host LF, so if the scatter in these relations were 
significantly smaller at z, it would increase r max . However, 
every observational indication, especially from the observa- 
tions in conflict with a low r ma v, is for an increasing scatter 
with z IShields et al l200l 120051: iTreu et alJl2004 IPeng et alJ 
120061 IWoo et al.ll2006l) . Moreover, if the mean M B H-host re- 
lation were higher in the past (z ~ 1), we still know observa- 
tionally that there is a non-negligible population (especially 
at the highest masses of interest) which has only passively 
evolved since such times (no star formation beyond the ~ 4% 
level, and Mrh can only increase) (see e.g .,|Kauffman n et alJ 
120031 IBell et dll2004aJ: iFaber et al]l2005l . so these systems 
must have been at or below the present relation at that red- 
shift (spheroid-spheroid "dry" mergers, which by definition 
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conserve the mean MeH-host mass/luminosity relations, can- 
not change this). Therefore, the scatter must increase with z 
in such a scenario, which gives a lower r max . 

Since the limits on T(z) rule out Mbh increasing strongly in 
the mean (i.e. for the general population) relative to the host 
bulge or total galaxy, it is next worth asking whether the dis- 
crepant T(z) from some of the observations in Figure[2]could 
represent some particular (short-lived) stage of development. 
The simplest possibility is that once an object has formed its 
bulge or spheroid, it will lie on the appropriate z = (passively 
evolved) relation, but that the observations see disks which 
have grown their BHs but have not yet grown a bulge. In such 
a scenario, T(z) in spheroids is given by the passive evolution 
tracks, and there is an independent, potentially much larger 
T(z) for "pre-spheroids." 

This pre-spheroid population, however, cannot constitute 
all disks or blue "star-forming" galaxies. First, at least 
some disks have survived since z ~ 1, and this would 
demand they have overmassive BHs relative to the total 
galaxy light/mass, where in fact the opposite is observed 
(i.e. Mrh correlates with bulg e, not total galaxy mass; e.g., 
Kormendv & Richstone 1995). Second, this scenario trivially 
modifies our calculation of r max - we integrate over spheroids 
with T(z) from the passive evolution tracks and determine the 
upper limit to T(z) in the remaining disk population. In Fig- 
ure |2] it is clear that the previous T max is already marginally 
consistent with the passive evolution tracks; adding an addi- 
tional source of high-mass BHs is not possible. Regardless, 
integrating down to a lower Lq (where disks begin to dominate 
the population) gives a similar constraint to those in Figure|2] 
since there is little room for "additional" BH mass density 
and, as noted above, the disk population evolves similarly to 
spheroids at the bright end (in the sense that the evolution is 
weak). 

The remaining possibility is that objects with BHs larger 
than our r max upper limits are an even more limited sub-class 
of galaxies. The natural expectation might be that AGN or 
bright quasars in particular obey this discrepant relation, since 
essentially all of the estimates of the high-redshift relation are 
from quasars. Locally at least, bright quasars constitute a tiny 
fraction of galaxies, so they can have much higher BH masses 
than their inactive counterparts (i.e. higher /iqso / n(z = 0)) 
without contributing much total BH mass. However, this 
phase would have to be short-lived, much like the quasar 
phase itself, or else it is equivalent to saying that a large frac- 
tion of non-spheroids obey this higher fi(z) and violates the 
integral constraints as above. Therefore, the host must "catch 
up" to have the appropriate stellar mass for the z = (passively 
evolved) relation in a short period of time. However, cosmo- 
logical cooling and infall generally requires a non-negligible 
fraction of a Hubble time, implying (if this were the mecha- 
nism to add bulge mass) that too large a fraction of galaxies 
would lie on the discrepant relation (given the integral con- 



straints above), and invoking subsequent mergers or interac- 
tions to transport gas would imply that the quasar somehow 
"knows" the system is going to interact before the interaction 
begins. 

This suggests that the mass required for the bulge to "catch 
up" is present, in the form of g as "about" to be turned into 
stars (see also e.g. Croton 2005), but this hypothesis encoun- 
ters a number of observational conflicts. To reach the appro- 
priate z = relation given the observations in Figure|2]which 
imply strong evolution in Mbh/M*, the systems would have 
to be > 75% gas at z > 1, which for the appropriate total 
mass range at these redshifts with these concentrations gives 
quiescent (lo wer- limit) star formation rates ~ 10 — 100's of 
M Q yr _1 (e.g.. lKennicutll998tlSpringel et al.l2005l) . The typ- 
ical quasar hosts would be late-type galaxies, and in the sub- 
sequent, brief period of rapid star formation the BH would 
be overmassive with the galaxy "catching up" (and could not 
significantly accrete without preventing this "catching up"). 

All of these predictions appear to be opposite the ob- 
servation s, with quasar spec t ra dom inated by posf-starburst 
activity ( Vanden Berk et al. 2006), typically early-type, 
bulge-dominated hosts (especially in the brightest quasars; 
Dunl onet alJ l2003t iFlovd et alJ |2004), and observed star- 
bursts usually showing high accretion rate BHs which are 
at that mo ment undermassive relative to their hosts (e.g., 
Borvs et al. 2005). Moreover, if star formation is ongoing in 
these objects to grow their bulges, it implies a lower mass- 
to-light ratio and suggests that the discrepancy between their 
MeH-host relations is even larger than that in Figure |2] am- 
plifying these difficulties. Furthermore, by the highest lumi- 
nosities and redshifts z ~ 1, the active fraction of BHs with 
masses M B h > 10 9 M Q is > 10% (relative to the tot al num- 
ber of such objects at z ~ 0; McLure & Dunlo p 2004), which 
creates some tension if these are a different population from 
the "finished" large spheroids at that redshift, which can al- 
ready account for the z = mass density of such large BHs 
(but is perfectly consistent if these are the same population, 
and quasars have large bulges). Finally, associating this dis- 
crepant T(z) with quasars connects it not with the cosmo- 
logical buildup of galaxies but with a particular evolutionary 
"snapshot" or event. This makes it difficult to explain why 
quasars at z > 1 would appear in bulge-less systems as a "sign- 
post" of bulge formation about to occur, while quasars of the 
same luminosity at z < 1 and low-luminosity AGN (which 
form a continuum in Lqso) would appear in normal, massive, 
bulge-dominated systems. 
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